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Abstract: We compute the value distributions of the eigenfunctions and spec- 
tral determinant of the Schrodinger operator on families of star graphs. The 
values of the spectral determinant are shown to have a Cauchy distribution with 
respect both to averages over bond lengths in the limit as the wavenumber tends 
to infinity and to averages over wavenumber when the bond lengths are fixed 
and not rationally related. This is in contrast to the spectral determinants of 
random matrices, for which the logarithm is known to satisfy a Gaussian limit 
distribution. The value distribution of the eigenfunctions also differs from the 
corresponding random matrix result. We argue that the value distributions of 
the spectral determinant and of the eigenfunctions should coincide with those of 



Seba-type billiards. 



1. Introduction 
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X. 

The study of quantum graphs as model systems for quantum chaos was initiated 
by Kottos and Smilansky JHj, |20| . who observed that the spectral statistics 
of fully-connected graphs are typical of those associated with generic classically 
chaotic systems. The relative simplicity of quantum graphs, together with the 
existence of an exact trace formula, has lead to the suggestion that their study 
might provide insights into some of the fundamental problems of quantum chaos 
|18j . This has motivated many works considering a variety of aspects of quantum 

graphs, El, 0, QH, B2I. EH, E2, ESI, ESI, EZI, EH- 

Studies |3] , 0] of a special class of graphs - the so-called "hydra" graphs or star 
graphs - have revealed spectral statistics that are not typically associated with 
quantum chaotic systems. These have been dubbed "intermediate statistics" in 
recent works, [HI, 0, CO] an d have been observed in a number of systems. We 
are motivated to investigate this further by studying the value distributions of 
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the eigenfunctions and the spectral determinant of the Schrodinger operator on 
quantum star graphs. 

A star graph consists of a single central vertex together with v outlying ver- 
tices each of which is connected only to the central vertex by a bond (figure ^) . 
Hence there are v bonds. We associate to each bond a length Lj,j = 1, .. .,v. 
We will often refer to the vector of bond lengths L := [Li, . . . , L v ). 



Fig. 1. A star graph with 5 bonds 



The Schrodinger operator on a star graph takes the form of the Laplacian 
— d 2 /dx 2 acting on the space of functions defined on the bonds of the graph 
that are twice-differentiable and satisfy the following matching conditions at the 
vertices: 

^(0)=^(0)=:<F, j,i = l,...,v 

= j* 

^(^■)=0, j = l,...,v. 

Here ipj is the component of the function defined on the j th bond of the graph, 
and ipj : [0, Lj] — > R with the convention that ipj(0) is the value of the function 
at the central vertex of the star graph. A is a parameter that allows us to vary 
the boundary conditions at the central vertex. 

The Schrodinger operator so-defined is self-adjoint, so there exists a discrete 
unbounded set of values < ko < fci < k-2 < • • • — ► oo such that fc 2 is an 
eigenvalue. It can be shown that fc = ±fc„ corresponds to an eigenvalue if and 
only if it is a solution of Z(k, L) = 0, where 

v 1 
Z(fc,L) := Vtanfcij - — . (1) 

K A 

We refer to Z(k,h) as the spectral determinant. Note that Z(k 1 h) has poles at 
k = (2n + l)ir/2Lj for each n G Z and j = I,..., v. The zeros and poles of 
Z(k, L) interlace. The usual definition of the spectral determinant would require 
these poles to be factored out. 

For simplicity we henceforth consider the case 1/A = 0. We shall employ the 
notation 

8Z 

Z'{k, L) = — (k, L) = V Lj sec 2 kLj. 

OK, 
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The eigenfunction corresponding to the n eigenvalue is found to be 

ip^fx) = A (n) C ° S kn ^ X ~ Li \ 
cos k n Li 



The constant is determined by the normalisation 

jr[ Lj \^\ X )\ 2 d X = i, (2) 

to be 



A {n) = 



y^yj— i Lj sec 2 k n Lj 



The value distribution of the eigenfunctions is determined by these normalisa- 
tion constants. For definiteness, we shall focus here on the maximum amplitude 
squared of the eigenfunctions on a single bond, 

Mn,L;v) := sup {\4 n) (x)\ 2 } 
xe[o,Li] 

= (A^secknL,) 2 

2 SeC fej-i-Z/j 



S/=i Lj sec 2 



We now state our main results. 
Theorem 1. For any fixed L > 

Hm ^ meas { L e [I, L + AL]<> : ^(Jfe, L) < = I £ ^d*, 

provided that kAL — > oo as fc — > oo. 

We emphasise that in theorem ^ we do not require that AL — > 0. In some later 
results we shall make this stipulation. 

Theorem 2. Suppose that the components of L are fixed and linearly indepen- 
dent over Q. Then 

If 1 1 1 r v 1 

lim — meas < fc € [0, if] : -Z(fc, L) < y ^ = - / -da:. 

K^oo if V J 7T J_ oc \ -\-x* 

Theorems ^ an d demonstrate the equivalence of taking a /c-average and a 
bond- length average at large k for the distribution of values taken by the function 
Z(k, L). Such a correspondence was noted in [2] for the spacing distribution of 
the eigenvalues of quantum graphs. 

In |17| it was shown that the value distributions of the real and imaginary 
parts of the logarithm of the characteristic polynomial of a random unitary 
matrix drawn from the circular ensembles of random matrix theory tend in- 
dependently to a Gaussian distribution in the limit as the matrix size tends 
to infinity, subject to appropriate normalisation. The distribution that appears 
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in theorems Q and [21 is known as the Cauchy distribution. It is related to the 
Gaussian distribution by the fact that both are examples of a larger class of dis- 
tributions known as stable distributions. Such distributions share the property 
that the sum of two random variables from a stable distribution is distributed 
like a random variable from the same distribution. Theorem ^ is a consequence 
of this fact. We also note that the density in theorems H an d|21 is independent of 
v when Z(k, L) is normalised as indicated. 

We next consider the distribution of values taken by Z'(k) when k = k n , 
n = l,2,.... 

Theorem 3. Let the components of L be linearly independent over Q. Then 
there exists a probability density P v (y), depending on L, such that 

lim l#(n 6 {1, ... , N} : \z\k n ,V) < R ] = [ P v {y)dy, 



with P v (y) = for y < 0. 



Theorem 4. For each v let the bond lengths Lj, j = 1, . . . ,v lie in the range 
[L, L + AL) and be linearly independent over Q. IfvAL — > as v — > oo then for 
any ReR, 

nR [-R 

/ P v (y)dy^ / P{y)dy 

J —oo J —oo 

as v — > oo. The limiting density is given by the continuous function 

' exp -V m(f)dC, y > 



where 



0, y<0 



i(0 : =^=e-« 2 / 4 + eerf(C/2). 

\/7T 



By comparison, the value distribution for the logarithm of the derivative of the 
characteristic polynomial of a matrix drawn from the CUE of random matrix 
theory, evaluated at an eigenvalue in the limit as matrix size tends to infinity, is 
Gaussian |T5) . 

The following results refer to the value distribution of Ai(n, L; v). 

Theorem 5. Assume the conditions of theorem^ are satisfied. Then there exists 
a probability density Q v (r]) such that 

lim e{l,...,N}: v 2 A t (n, L; v) < R} = [ Q v (v)dr) 

JV^oo iV J a 

where the density Q v {rj) is independent of the choice of bond i but depends on 
L. 
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Theorem 6. Assume the conditions of theorem^ are satisfied. Then for each 
R>0 

Q v (r))dr)^> / Q(r])dn. 
10 Jo 

as v — > oo. The limiting density is given by the function 

which is continuous on (0, oo). Herem(t;) is as in theorem^ Q(i]) has asymptotic 
expansion 

V2 r e-« 2 / 4 B/a . 



as r\ — > oo. 



VLir 2 n 3 / 2 J-oc m{€) 



The proofs of theorems rely on an equidistribution result of Barra and 
Gaspard j^j- We review this work in section |3J 

The limit v — > oo is analogous to the semiclassical limit h — > [20] . We note 
that theorem Ej describes the wave functions on a vanishingly small fraction of 
the graph. It thus goes beyond the information provided by the Schnirelman 
theorem. It instead corresponds to the Gaussian value distribution for the wave 
functions of classically chaotic systems implied by the random wave model 

The value distribution of the eigenvector components of asymptotically large 
random matrices are particular cases of the \^ density 

where the parameter (3 takes the values 1, 2 and 4 in, respectively, the orthogonal, 
unitary and symplcctic ensembles (see for example JHl)- When 0=1 the density 
is called the Porter-Thomas density. It is characterised by 0(?7 -1 / 2 ) behaviour as 
7] — > and 0(r/ _1 / 2 e - ''/ 2 ) as t] — > oo. The limiting distribution we find in theorem 
El completely determines the value distribution of the star graph eigenfunctions 
(see the appendix) and has a significantly different shape (c.f. equation J3J and 
figured below). Other quantum systems for which the value distribution of the 
eigenfunctions has a non-random-matrix limit are the Cat Maps 22 . 

In [3] a correspondence was noted between the two-point spectral correlation 
functions for star graphs and a class of systems known as Seba billiards. The 
original Seba billiard |28) was a rectangular quantum billiard perturbed by a 
point singularity. More generally, we describe any integrable system perturbed 
in such a way as belonging to the same class [25] . We conjecture that the results 
derived in the present work will also apply to systems in the Seba class. 

The remainder of this paper is structured as follows. In section [21 we prove 
theorems ^ an d El I n section we treat the finite v cases, theorem [3] and [5| 
In sections 01 and we prove, respectively, theorems 0] and El developing the 
necessary machinery in sectional SectionEJis devoted to numerical computations 
that illustrate our results. We develop more fully the connections between the 
present work and Seba billiards in section [7| 
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2. The value distribution of Z(L, k) 
Lemma 1. Let L > and £ be real constants, then 

Y pL+AL 

lim — — / exp(i£ tan kL)dL = e~ 

k^oo AL J I 

uniformly for kAL — > oo as k — ► oo. 

Proof. By the periodicity of the integrand we may shift the range of integration 
by multiples of n/k so that without loss of generality we may take L in the range 
< L < n/k. We write AL = n(n+p)/k where n £ Z and < p < 1. Then by 
the periodicity of the integrand, 

i>L-\-izn / k+irp / k 

I exp(i£ tan kL)dL 

pL+pir/k r^/k 

= / exp(i£ tan kL)dL + n / exp(i£ tan kL)dL 

Jl Jo 

r-L+p-K/k n rCC e i£ z 

= / exp(i^tanfci)dL + — / ^dz 



where the substitution z — tanfcL has been made. We note now that n/kAL 
7r~ as k — > oo and 



^ i-L+pTv/k 

-^jr j exp(i£tan/cL)dL 



7T 

< — ; > as k 

~ kAL 



A simple application of Cauchy's residue theorem allows us to evaluate the in- 
tegral 

dz = 7rc- |cl . (5) 



oo e i(z 



-oo 1 + Z 2 

□ 

Proof of theorem^ We use here the characteristic function. With bond lengths 
chosen from a uniform distribution, 

j rL+AL r-L + AL 

E L (exp(iCZ(fc,L))) = jALrJ_ -"J exp(iCE; = itanfcL,jdL 1 ---dL.„ 

/ 1 f L+AL \ v 

= I -^j? J exp(i£tanfcL)di j . 

The subscript L indicates that the expectation is with respect to an average over 
bond lengths. Since the map t t" is continuous, lemma ^ together with JSJ 
allows us to deduce that 

lim E L (exp(iCZ(fc,L))) = e^. (6) 

k — MX) 



Eigenfunctions and spectral determinants of star graphs 



7 



The limiting density corresponding to the characteristic function on the right 
hand side is given by 



The theorem follows now from the classical continuity theorem for characteristic 



The proof of theorem|21uses Weyl's Equidistribution theorem. This celebrated 
result has numerous applications in analysis and number theory. We state 
here the form most convenient for application to our current work. Let T" be 
the i>-dimensional torus, with sides of length tt. 

Theorem 7. Let f S C(T V ), and let the components ofL be linearly independent 
over Q. Then 



where dx = dxi ■ ■ ■ dx v denotes Lebesgue measure. 

We shall use Weyl's theorem as our main tool to relate fc-averages to bond 
length averages. It is for this reason that it is crucial that the bond lengths are 
incommensurate . 

We remark that theorem can also apply to more general functions such as 
piecewise continuous functions through an argument similar to the one in the 
following lemma. 

Lemma 2. Theorem^can also be applied to the function 



Proof. We treat the real and imaginary parts of / separately. The functions 



are smooth everywhere apart from at an essential singularity when x% = it/2 
for some i, which we tame in the following way. Let e > . We can construct 
functions <fr and ip satisfying the conditions of theorem [7| such that 




functions ( |14| chapter XV). 



□ 







(7) 



This implies 
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From (0 and theorem 

-1 / ^(x)dx<kminfi [ h(kL)dk 

< lim sup — / fx(kL)dk<^- [ 0(x)dx. 

K^oc & Jo K v J r v 

The ends of this inequality differ by e which can be made arbitrarily small, so 
we see that hniA'^oo K~ l J K /i(fcL)dfc exists and is equal to 



/i(x)dx. 



1 

The extension to ji and hence / is obvious. □ 



Proof of theorem We begin in the same way as in the proof of theorem ^ 
In this case k is chosen uniformly from the interval [0, K] with K > 0, so that 
the characteristic function with respect to this uniform distribution is 



E K (exp(i(Z(k,L))) = J exp (iC SjL x tan kL j) dfc - 
By lemma [21 we can write this integral as an average over the torus as K — > oo : 
lim E K (exp(iC^(fc,L))) = — / exp (iCELi tana:,,) dx 

= — / cxp(i£tanx)da; 
V 71 " Jo J 

Following the substitution z — tana; in this final integral, we have 



lim E K (exp(i(Z(k,L)) = c 

K—>oo 



and the theorem follows from the same arguments used in the end of the proof 
of theorem ^ D 



3. An equidistribution theorem 

Barra and Gaspard [2] observed that the condition for k to be an eigenvalue of 
a graph can be written in the form 

G(fcL) = 0, 

where G is a function that is periodic in each variable. For star graphs G is 
defined on T v by 

G(x) = taniri + • • • + tanx v . 

The equation 

G(x) = 
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defines a surface S embedded in 

fc (x o ) = 



y. A flow can be defined on TP by 

xq + kh (mod 7r). (8) 



Since k = is an eigenvalue for star graphs with Neumann boundary condi- 
tions considered here, we take Xo = in this case. 

At each value k — k n we have an intersection of this flow with the surface S. 
We note that the angle between the normal to the surface £ and the flow cf> k is 
given by 



|L- VG\ 



cos ft = 



• SGC Xd 



INIIvgii 

For star graphs, 

VG(x) = (sec 2 a;i, 

Hence there exists a constant ci > such that cos 9 > c\. This means that the 
angle between the flow and the surface £ is uniformly bounded away from 0. We 
can therefore parameterise S locally by v — 1 real variables £ = (£x, ••• ,£,v-i) 
so that for x G S, 

Xi = Si (£) 

and G(si (0,... ,«„(€)) = 0. 

The central result of Barra and Gaspard is the existence of an invariant 
measure on the surface S. 



Theorem 8. Let f be a piecewise continuous function S 

N 



Then 



where the measure v is given by 

di/(e) 



k mm 

and d£ = d£i • • • d£„-i is Lebesgue measure. J is the Jacobian determinant 



(9) 



(10) 





L v 


ds± 






d£i 


dsi 


ds v 




dtv- 



For completeness, we sketch a proof of theorem|Hlfor star graphs with v bonds. 



Proof. Let / : TP 



be an extension of / to T, so that / 
/(x) = /(x) for all xeS. 



/, i.e. 



We let / be constructed in such a way that for all £ S S, f{4> k {£)) is a differen- 
tiablc function of k with compact support in some neighbourhood of k = 0. 

Let e > 0. We construct the set which is a thickening of S in the direction 
of the flow 6 k , 



£ e>L := {x 6T":](e£i£ [-e, e] : x = fe (O} C TP. 
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We define 1^, the indicator function of a set A, by 

M*) ■■={llt A A 

The indicator function TL^, L ( x ) * s piecewise constant. 

By the differentiability properties of / we can write for every x££ 

/(x) = /(x) = i | £ /(fcL + x)dfc + 0(e) (11) 

as e — > 0. The implied constant does not depend on x. Setting x = fc„L gives 

N n e 

- ^ /(fe n L) = — ^ / /((fc + fc„)L)dfc + 0(e). 

n=l 6 n=l^ _e 

Let the mean density of zeros of Z(k, L) be d: 

J:= lim (12) 

Then 

J im v E /( fc » L ) = Tl J im 1? / Z(fcL)lU, L (fcL)dfc + O(e) 
N^oo Jy — ' Zed k^oo K In 

n—1 

= 4l[ /(x)li; e . L (x)dx + 0(e), 

applying theoremdto the piecewise continuous function /(x)ls t L (x). Changing 
to the system of coordinates (t, £) on T" via the change of variables 

%i = Lit + Sj(£), 

gives 



n=l 

Since this is true for all e > 0, we deduce that 

iV 

i 

im — 

JV 



1 w . . 

} im at E /( fc « L ) = \ /(^) J (0d€- (13) 



By setting / = 1, we see that 

d=\ I J(€)d€, (14) 
to complete the proof. □ 
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We note incidentally that d can be evaluated using spectral methods j2D] to 
give 

1 v L 
d = — L.j = - h O(vAL) asvAL^O. 

7T ^ J IT 

3=1 

We observe that the right hand side of equation © can formally be written 
in the form 

^ /(€)di/(€) - J jf • • • Jf /(x) [L • VG(x)]e i « G dx 1 • • • d^dC, (15) 

where now / is a function TP' — * R of an appropriate class. This follows from 
writing J„ /(£)dz/(£) in the equivalent form, 



(16) 



(17) 



~dj s j /(<,0[L ■ VG]5(G(t,0)J(C)dM^. 
We then write the (S-function as the limit of the sequence 



S m (x) 



1 

2^ 



e iCa; dC 



sm mx 



as m — > oo, and changing back to the usual Cartesian coordinates on T v . We 
need to show that 8 m is an appropriate <5-sequence for the function / that we 
consider. We will check this point directly in the calculations where the identity 
(|15|l is used. We shall also justify taking the ^-integral in l|17|) outside the integral 
over T v . 

Proposition 1. For a star graph with v bonds and the parameterisation 

Si=£i i = l,...,v-l 
s v = — tan -1 (tan£i + • • • + tan£„-i) 

J(£) takes the following form 

Li sec 2 £i H + L v _x sec 2 £„_! 



Jit) = 

l + (tan^i + --- + tan^_i) 2 
Proof. Differentiating gives 



(18) 



dsi 



for i < v and 



ds v 



sec 2 



d£j l + (tan^i + --- + tan^_i) 2 ' 
For ease of notation we write D := 1 + (tan£i + • • • + tan^„_i) 2 . Thus we have 

L v —i L v 

q - sec 2 g 







Li 


L 2 ... 




1 


••• 


m = 





1 ■•• 







••• 



— sec g 2 



-sec 
D 
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To complete the proof we employ the identity 



Otl 


a 2 ■ 


• Otn-i 


«n 


1 


• 


■ 


ft 





1 • 


• 


02 





• 


• 1 





= {-IT ^-«n + £«*^ » 



which may be readily checked by induction. 



□ 



We note in passing that the explicit form of J(£) given above together with 
theorem[H]provides a convenient representation for use numerical studies of eigen- 
values and eigenfunctions because the zeros of Z(k, L) do not need to be com- 
puted explicitly. 

Proof of theorem We take as the function / in theorem [H] 



/(x) = 1 



(-00, fl] 



3 = 1 



Then we define P v (y) by 



tan£„_i)) 



Pv(y)dy = — 1 •••/ /(£i,...,& J _i,-tan~ 1 (tan£i + - 

-oo n a JO JO 

xJK)dfi-dJ,-i, 

where J(£) is defined by l|18|l . Since sec 2 x > for all x € R, it follows that 
P v (y) = for y<0. □ 

Proof of theorem In this case, we take as the function / in theorem |H1 

/ 2v 2 sec 2 Xi \ 



/(x) = 1 



[0,R] 



Ej . sec2 ^ 



We take as the parameterisation of S 

s j =€j> l <3<i 
Si — — tan (tan £i + 
s j = £j-i> i <3 <v. 



tan£„_i) 



This does not change the form of J(£) from that in JTSJ, but introduces extra 
symmetry in /. Then define, as before, 



Since /(s(£)) is symmetric in £i, . . . ,£ v -i we see that P v (rj) is independent of 
the choice of bond i. □ 
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4. Value distribution of Z'(k n ) in the limit v — > oo 

We take as the function / in (|15|l the characteristic function for the value dis- 
tribution of the derivative of Z(k,TJ), 



Z'(k,L) =J2 L i se ° 2 kL j- ( 19 ) 
4=1 

Since we are stipulating that vAL — > as v — > oo, we replace Lj by X where it 
does not multiply k in 119fl and take as our / 



/(x) = exp ( ^ sec 2 



4 



Let the quantity in which we are interested be denoted E v {(3). Then 



£„(/?) := lun^ 1 ^ exp ( £ sec 2 fc^ 

n— 1 



2irdvJ_ 0O ir* J io \4" 



X 



exp I — — ^- sec 2 £j + — tanxj J dxd£. (20) 
4=1 V U ^ 4 



We have made the re-scaling £ i— > C,/v. This is a natural normalisation since 
Z{k, L) is a sum of w terms. 

We exploit the symmetry in the integral in (|20|l to write 

1 f°° 

where we have replaced d by Lv/tt and defined the integrals 

Ji(/3, £):= — / sec 2 a; exp [ — sec 2 x + — tana: J da: 
Ti" Jo V w v J 

and 

^2 (A C) : = ~" / exp (—-^- sec 2 a; + — tana; ) da;. 
■k Jo \ v z v J 

We note that I2 is uniformly convergent in £ but that Ji is not. 
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The integrals I\ and li- The substitution z — tanx gives 

MAC) = - r exp f-^(l + z 2 ) + dz 
v 1 / C 2 



quoting a standard integral. 

That 6 m defined in (|17f) is a (5-sequence for the function exp(iaz 2 ) follows 
from the equality 



1 

2^ 



CO /"CO 

e iaz2+ic( - z - w) dzd( = expfW 2 ), 

oo J — OO 



which can be checked by direct evaluation of the integrals. This and uniform 
convergence of I 2 in ( justifies the operations that lead to identity (|15|l . 
We can treat 7 2 in a similar manner to that in which we treated I\ . 

7T J_ co \ V 2 V J 1 + Z 2 

We first write 

1 1/1 1 



1 + z 2 2i \z — i z + \ / 
so that J 2 can be decomposed into a difference of two similar integrals 

/ 2 (/3,C)=:/ 2 -(/U)-/ 2 + 03,O- (22) 

Observing that 



exp(- — (l + z) + —)=exp [~—~ m ~— [*+ m 



we can write 

1 f°° ( [ P L ,i 2^ '< z \ dz 
h(0,O-=ir-J exp -^(l + z 2 ) + ^ ; 

27T1 J_ oc \ V 2 V ) Z — I 

1 ( ifiL iC 2 \ f°° eM-^Ly 2 /v 2 ) 



~ 2th CXP v 2 + 4(3L ) y - i - (v/2(3L dy 

via y = z + C,v/2(3L. We make the change of variable 

r = y. 

v 

This is permitted, because it rotates the contour of integration into the second 
and fourth quadrants of the complex plane, where the analytic function e~ lz 
decays rapidly. 
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In the case (v/2(3L > — 1 the new contour of integration avoids the pole at 
y = (v/2f3L + i (figure EJ) and Cauchy's Theorem yields 



h (A = ex P 



iJ3L , iC 2 ^ 

( ,2 



-dr. 



-^(i + ^1 

« v ^ 2/3L'' 



This integral is standard, and may be found in, for example, [1, (Equation 7.1.4) 

-di = e _z erfc(— iz), for 3m z > 0. 



z - i 



The result we get is 



h W> 0=2 ex P 
x crfc 



1 



1 / ipL iC 2 \ / i/3£ 
= — exn I - — I =■ 1 exr) I — 



' 4/3L/ 



1 



exp 

i( 



2/3L 



= — exp — erfc 



2(3L 
C 



2Vi/3L 




If £v/2(3L < — 1 then the contour encloses a pole (figure |3J). In this case, 



cxp(— \j3Ly 2 /v 2 ) 
y-i- Qv/2I3L 



dy = 2niR 



-dr 



'ij3L ,. 
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where R is the residue at the pole 



R = exp 



\f3Ly z 



y=i+Cv/2pl 



iC 2 ^ iPL ( . 

A(3L ir v 



so that we also get in this case 



I2 (0,0 = ^cxp( i lerfc 



2\/i/3L 



Treating in a similar way, yields an expression for I 2 , 



/ 2 (/U) = ^/"eric 



C + Vi^) + I e -c/, erfc 

2y/i/3L v I 2 



< 



2Ji(3L 



+ ■ 



(23) 



4.2. Some estimates. 

Lemma 3. Let —\Jv < ( < ^/v, then 

wow,or'-^»p(^)*»o 

x (l + 0(t,-" 2 ) + 0(C 2 /f)) (24) 
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as v — > oo, where 

M(PX) :=exp 



2-v/i 



i/0£ / -C 2 \ 
= — exp I =■ I 



Cerf 



c 



2 x /i/3L / 



Proof. A key step will be to make a uniform expansion of i2(/3, C) as 
By Taylor's theorem, 



c 



c 



exp ±- = 1 ± - + 0(C u ) as u -> oo 



A second application of Taylor's theorem yields 



erfc 



7 i/3L 



erfc 



2 Vif3L 



2sJ\j3L v J 
where the remainder term is 



■ exp 



4i/3L 



R 1 (0 = -i(3L I — erfc(z) 



L & 
d^ 2 



(1 - er)der. 



.+<T- 



The second derivative of erfc is 

d 2 



This allows us to estimate 



\Ri\< 



< 



—— erfc(z) = —=ze~ 
dz z Wtt 



ICI 



A/3L 

Jo \2y//3L 

m lcl/v ( JO. 



+ o( v - 1 ) 



V 7 !" \2VI3L 
= 0(v / w) as w — > oo, 
uniformly for |£| < v 7 ^- Substituting JSHJl and (|2"?jl into (|2"3)> gives 



wc)=i-^#e XP r^;v^rf 



c 



\4i/3Ly « \2y/ii3LJ 
where i?2 is a combination of the errors in (|26|) and So 
|i? 2 (C)| =0(1 + C 2 ) asw^oo. 

We thus have that 



= — exp I I 



\4i0Lj 



£ erf 



c 



MO 



2Vi/3L 



< 



2v^I , |C| , \R 2 \ 



7TV y/v W 3 / 2 

< 2 for all ( 6 [— y^, ar± d u sufficiently large 
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We are now in a position to understand the asymptotics of {Ii(f3, Q) v 1 • We 
first consider 



In 



(32) 



where 



1^3 (a) I < 



l-\a\/v 



provided that \a\ < v. We take 



c 



MO 



2Vi/3L 



which satisfies |a|/ < 2 for v sufficiently large by Ij31|l and a = 0(|C|) as 
|C| -> oo for |C| < b y (EOJ)- We see that 



\R 3 (a) 



0(\a\ 2 /v) as v — > oo, 



so exponentiation of lj32J) gives 

(l - J) ^ = e- (l + O^ 1 / 2 )) (1 + 0((1 + C 2 )/-)) 
where the error estimates are uniform for |C| < y/v as u — » oo and 
2 x /i/3Z /-C 2 \ . „/ C 



(33) 



exp 



exp — c erl 



We can also write 

H/3,0 



v 1 



: exp 



I < 2 \ 



[f + o((i + c 2 )/«)] 



(34) 



Thus, taking 113*4")) together with l2l"}) and (f""""j) gives the required estimate. □ 



Lemma 4. Let \C\ > y/v, then 

|7 1 C9 > C)(/ a C9,C))"- 1 l <^ + +OK" 3 ) 

Proof. We make use of the asymptotic expression 

erfc(z) = ^e~ z2 (- + 0(z" 3 ' 
as z — > oo, valid for | argz| < 3tt/A. We see that 



as |d — ► oo. 
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Similarly, 



e- c/v erfc 



-c 



1 



2Vi/3L 



< | y/i0L 



2Vi/3i 



o(c 



Adding these gives the estimate 

i«ai<f(e4 

Taking together with the estimates | J2 C) I < 1 an d 

l^i (P, 01 < 



+ 0(C 3 ) as ICI ^00. (35) 



gives us the estimate we require. 



□ 



Proposition 2. With the notation above. 



lim — 

v— *oo 2l? 



/!(/?, c)(/ 2 (/?, or -1 dc 



1 



x 1 f -c 2 



: exp 



Proof. We split the region of integration as follows 



-/ / 1 (/3,C)(/ 2 ( / 3,C)r- 1 dC=- / h(J3,Q{h(P,Q) 
2w J-00 2u J_ vT j 



M(/3,C)dC. 



2i; 



The function Af is bounded and satisfies 



ICI>v^ 



JiC/3, C)(/ 2 (/3,C))"- 1 dC. (36) 



|M(/3,C)I <exp 



2v^i ICI 



for ( sufficiently large. Integrating (|2~lf> gives 



2 v 



f -c 2 



M(/?,C)dC 



O 



(37) 



and the integral in the remainder term converges as we let v — > 00. 
To deal with the integral 



1 



\(\>Vt 



/iCAcxwor-'dc 
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we make use of the result of lemma 0J giving 



v J\c\>vt 



h(J3, C)(/ 2 (/3,C))^ 1 dC 



< 



1 



0l + e 



as v — > oo. 



d( + 0(iT 



(38) 



Hence substituting (|37|l and l|38|l into (|36|l and taking the limit v — > oo gives 
the required result. □ 



^.3. Properties of M(/3, £). We wish to rotate the variable £. However, we need 
to check that the function M does not blow up for large |£|. 
If C = Re w , then 



erf 



erf 



R 



2y/ (3L 



\2y/i(3L J 

as R — > oo, provided that < < 7r/2 p. So, 



3 i(9-V4) =1 + 0(7?-!) 



(39) 



exp — £erf 



2\/i/3L 



= exp ^-i? SHc 

— J? COS 0, 



e 1 " erf 



C 



\2 v /i/3L / 

= e - — J 0(l) 

as R — > oo provided < 8 < tt/2 



and convergence is exponentially fast. Similarly, 



/ 2Ji/3L 
exp — — exp 



4i/3L 



{ -2V(3L 
exp — — exp 



R 2 



A/3L 



■ sin 26 



I R 2 nn * 
x cos — =■ cos 20 H — 

\AI3L 4 
1 as R — > oo 



provided < < tt/2. 

Hence, if < 6 < it/ 4 then 



lim \RM(f3,Re ie )\ = 



and we can make the change of variables £ = £yif3L: 
f°° 1 



/ 4 M(/3,^i/3L)de (40) 
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We observe that 



M[-i.i,\/llL) <-xp( — ^=^8Ze-« 2/4 - £^/i/?£erf ( | 



exp — yi/3Lm(£) 



where 



i(0 :=^c-« 2 / 4 + eerf(^/2). 



Proof of theorem^ m satisfies the bound m(£) > 2/y / 7r, so M(/3, £y/i/3L) 
is bounded for all /3. By the Weierstrass M-test the integral in (|40|l is uniformly 
convergent and hence 

lim E v (l3) 

v — >oo 

is a continuous function of (3. We appeal, once again, to the continuity theorem 
for characteristic functions to deduce that the limiting density P(y) exists and 
is given by 

"1 7*00 /*QQ I 

P W = TTJ2 / e-« 2 / 4 M(/3,eVi/3L)e^dCd/3 

J cc J CO 

i /"CO /"OO y 

= 2^ me J J e ~ e/AM iP^y^ L )^^- ( 41 ) 

The integrand is dominated by 



exp \-?/4 - sJttL[3 

so Fubini's theorem allows us to switch the order of integration. We quote the 
standard integral 

1 e a* + ^ dx = I _ A A /Z exp ( Z^) erfc ( - b 







a 2a V —a \ 4a / V 2J —a 



valid for [He a < and use this to perform the (3 integral in (|41J) . This leads to 
the result 

^) = J** (fT - -f^J m(e)6rfc l^^iT- J d? ' 

which reduces to the form given in the statement of the theorem upon noticing 
that SHe{erfc(i<9)} = 1 for all 9 € M. □ 
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5. Value distribution of the eigenfunctions in the limit v — > oo 

To prove theorem [S] we use a standard approximation argument. We introduce 
the smoothed (S-function 

S e (x) := - / e-^l +i ^d/3 = (42) 

Proposition 3. Let Q v (rj) be related to Q v {rj) by 

Qv(v) : = \Qv ( - 
For any fixed e, 

poo poo 

lim / S e (ri-rf)Q v (Tf)dif= / 8 e (r, - rf)Q(rf)drf, (43) 
v ^°° Jo Jo 

uniformly for rj in compact intervals, where 
Proof. Let 

AT 



f°° 1/1 

Q e , fl (»7) := / « e fa " tfQvWW = lim _ £ J e U - 

u n— 1 x v 



n, L; u) 



We use the identity 

<5 £ ^-^) B,\ Li u.h r M 

with 



A = — ^ sec 2 kLj 



and 
Thus 



v 

3=1 



B = 2 sec 2 fcLj 



p / ^.Lsec^kLA 2sec 2 fcL 4 f°° f f L . ., 

* e U - =^ TTT^ = / exp i/3 277 - — + 2ie sec 2 fci, 

\ sec^ fcLi / 7T 7n V V v 



x J^Jcxp ^— sec 2 kLj^j d/3. 



Applying identity Ijl5|l to this function and following the method in sectional 
leads to 

c\ j r-OO pOO 



ir 2 dv 



o J -oo 

+ (v- l)/ 3 (/3, C)^i(/3, C)(/ 2 (/3, C))"" 2 dCd/3 =: Afo) 
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where the new integrals are 



h(P, C) : = — / sec 2 a; exp ( i/3 ( 2rj — 4j + 2ie ) sec 2 a; + — tana; ) da; 



and 



sec 4 a; exp I i/3 I 2r] t 2k I sec 2 x + — tan.c ! dr. 



J3 and Z4 converge uniformly in £ and r\. Integral J 3 closely resembles integral 

h. 



hW,0 



By making the substitution z = tana;, I4 reduces to 



+ O e (v 1 ) as v — ► 00. 



(44) 



1 



iCz 



J 4 (/3, C) = - J (1 + z ) exp I ij9 (27? - — + 2ieJ (1 + z z ) + ^ dz 
— O e (l) as u — > 00. 

This estimate ensures that Q e ,v(v) is dominated by the second term. Since I3 is 
also bounded as v — > 00, the analysis of proposition [3] holds and 



lim QcvM 



v— >oo 

2 



1 



OO /"GO 



7T 2 V2e - 2177 J J^oo /3ViL 
1 



1 / c 2 
exp 2i/3?7 



4i/3L 



2e/3 M(/?,C)dCd/3 



4^e 

-7T-Z 



OO /"OO 



, , exp (-L- JipLmiO + 2i/3,y - 2ep] ^§df 

V2e - 2iT7 Jo V 4 V / VP . 



where we have made, once again, the substitution £ = ^^J\j3L in the final line. 
Putting 2 = (3 reduces the /3 integral to 



d6 



■. exp ( - ^ m ^\ ) ei-fc 



2^2e - 2i?7 \ S^ + ie ) 



Lm(0 



2V2r? + 2ie / ' 



lim Qe,v(ri) 

v — >oo 

1 



2tt 3 / 2 



3m 



1 



77 + ie 



exp 



e ^) 2 V-,-f y : 

4 8(77 + ie)7 I 2iV2?7 + 2ie i 



5 e (r? - r]')Q{rf)drf 



□ 
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We note that Q(tj) is a continuous probability density on (0,oo). 
Proof of theorem® Let < a < b be fixed. Then let 

^■.qfo) : = / S £ (r]-y)dy. 

J a 

lj Q b ](r]) converges pointwise as e — * to the function U[ a everywhere apart 
from at the end-points a and 6. Given cr > 0, consider the function 

Xl(v) ■= l[ a -a,6+ CT ]( J ?)+ CT ' 

There exists e > such that:- 

— < Xiiv) — 2c for ?7 < a — 2cr and 77 > 6 + 2<r, 

— 1 < Xiiv) < 1 + 17 f° r a < ?7 < 

— < Xi (^) < 1 + ° for a - 2cr < r) < a and b < r] < b + 2a. 
This construction is illustrated in figure 01 




Similarly, the function 

X2{ri) = l[a+ CT ,b- CT ]( ? 7) -° 

satisfies for e sufficiently small:- 

— — cr < X2(ji) < for rj < a and rj>b, 

— 1 - 2cr < X2{rj) < 1 for a + 2a < rj < b - 2cr, 

— —a < X2(v) < 1 for all a < 77 < a + 2cr and b — 2a < r\ <b. 



Eigenfunctions and spectral determinants of star graphs 



25 



So that for all r\ S [0, oo) 

Xa(v) < 1 [«,6]( ? 7) <Xi(v)- (45) 

Also, 

/>oo />a+2er 

[Xi(v) - X2(V)] Q(v)dv < 3fr / Q(ri)dr] + (1 + 2a) / Q(r ? )dr ? 

Jo J a-2o 

+ (1 + 2a) f + Q( V )d V: 

which can be made arbitrarily small because Q is a continuous probability den- 
sity. It follows from proposition [3] that 

/>oo />oo 

lim / Xi(v)Qv(r])dr] = / Xi( 7 7)Q( 7 7)d?7 
t, ^°° Jo Jo 

and similarly for \2 ■ Hence, we can use the argument of lemma |21 mutatis mu- 
tandis, to deduce that 

lim / Q v (v)drj= [ Q(v)d V . (46) 

Making the substitution n h- ► \jr\ then completes the proof of convergence. 
Expanding the error function in @ as 



erfc ( ^*f>) = -1= exp ( ( ^ + Ofr"*) 



where the implied constant does not depend on £, yields 



^fa) = ^2 + 5/2 ) as ry -> oo, 



where the constant b is 



V2 [°° e-« 2 / 4 ^ 0.348 



VZtt 2 J-oo VI 



□ 



The algebraic decay of Q{rf) is in contrast to the exponential decay of the \\ 
density. 
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6. Numerical Results 

The results presented above show close agreement with numerical computations. 
We present these computations now by way of illustration. 

In all the figures in this section, the choice of L = 2 has been made. 

Figure El shows a comparison between a numerical evaluation of values taken 
by the spectral determinant and the Cauchy distribution. The numerical eval- 
uation was based on a star graph with 7 randomly chosen bond lengths, and 
100,000 samples of k. 

Figure El shows a comparison between the distribution of values taken by the 
derivative of the spectral determinant at its zeros, and the corresponding nu- 
merical evaluation. Plotted is numerical data for a 70-bond star graph, together 
with the v — > oo limiting density given in theorem 01 Once again we see good 
agreement. 

In figure we compare a numerical evaluation of the density of values taken 
by the maximum norm of eigenvectors of a 50-bond graph to the v — > oo limiting 
density given in theorem Also plotted for comparison is the density of the Xi 
distribution associated with the COE of random matrices. 



7. Connections with the Seba Billiard 

The correspondence between the spectral statistics of quantum star graphs and 
those of Seba billiards with periodic boundary conditions has already been noted 
3 . This is due to the fact that the spectral determinant for the star graphs 
may be re- written in a form similar to the spectral determinant of a Seba billiard: 



OO 

Z S cba(£) = -75T 



where the are the energy levels of the unperturbed system. Both spectral 
determinants have infinitely many poles of first order, which separate the energy 
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Numerical simulation - 
Predicted density - 




Fig. 6. The value distribution of Z'(k) 



Numerical simulation - 
Prediction - 

COE ■ 




Fig. 7. The value distribution of Ai(n, L;t>) 



levels of the perturbed system. We therefore expect the value distribution of 
the spectral determinant of a Seba billiard to be Cauchy to be consistent with 
theorems ^ and |21 

This conjecture is supported by figure|Hlwhich is a plot of the density of values 
given taken by the function 



K 



ir(d) 



(47) 



for K = 3000 unperturbed levels of a rectangular quantum billiard with Neu- 
mann boundary conditions, with E distributed uniformly between E^ m and 
i?2oo - The constant (d) is the mean density of levels of the system and it takes 
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the place of the constant v 1 in theorems Hand El The fit to a Cauchy density 
is convincing. 



Numerical simulation - 
Cauchy distribution - 




Hi™ 



Fig. 8. The value distribution for the spectral determinant of a Seba billiard 



If we treat the unperturbed levels in 14711 as independent identically dis- 
tributed random variables with a uniform density then the random variables 



E. 



(o) 



E 



have a distribution that falls into the domain of attraction of the stable Cauchy 
density. That the limiting density is Cauchy is then a classical result of proba- 
bility theory [TH . 

We now present an argument which suggests that the normalisation constant 
associated with the wave functions of the Seba billiard also shares significant 
features with the normalisation constant of the star graphs J2J. 

The wave functions of a general Seba billiard can be written in the form 



V>„(x) = A, 



E 

k=l 



4 0) (X )4 0) (X) 



(48) 



where E n is the n th energy level and Ejp and tp^ are, respectively, the energy 
levels and wave functions of the original integrable system of which the Seba 
problem is a perturbation. The Berry- Tabor Conjecture 7, asserts that the un- 
perturbed levels Ej?' are distributed like Poisson variables; that is independent 
and random. We fix the usual normalisation 



|^ n (x)| 2 dx= 1, 
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which leads to a value for the constant A 



-n 




Neumann boundary conditions and sides of length a 1 / 4 and a x / 4 the wavefunc- 
tions are 



V'i°L(^l/) = 2cos(^)cos(m^a 1 / 4 ), n, m = 0, 1, 2, . . .. (50) 



If we position the scatterer at the origin, then |V4,m(x )| = 2. This billiard 
problem is equivalent to the billiard with periodic boundary conditions desym- 
metrised to remove degeneracies in the spectrum. Provided that the constant a 
satisfies certain diophantine conditions (see |23U13| for details) then A\ in H49JI 
is the reciprocal of a sum of functions with poles of second order distributed 
independently. These poles play the role of the singularities of the functions 
sec kLj which appear in the normalisation of the quantum graphs. Such poles 
determine the rate of decay of the tails of the relevant probability distributions, 
and this implies that the analysis performed in the present work also holds for 
this billiard problem. In particular, we conjecture that the distribution of the 
square of the i th coefficient of the eigenfunctions in the basis ) is the same 
as the limiting distribution of Ai(n, L; v). 

We present in figure the distribution of values taken by 



where n is now a random variable uniformly distributed on {1000, . . . , 2000} and 
we take K = 3000, i = 1500 and a = (s/5— 1) /2. The constant c, which in general 
may be expected to depend on K and the distribution of n, is required to ensure 
that the sum of terms in the denominator is normalised and to compensate 
for the fact that the functions are not periodic. In order to compare with the 
corresponding results for star graphs we require that the tail of the distribution 
of c(E^ — E n )~ 2 is asymptotic to the tail of the distribution of (2/Z) sec 2 k n Li. 
Assuming n to be distributed between n max and n m ; n , a heuristic examination 
of these densities leads to the association 



where £ max and E m i n are respectively the energy levels corresponding to n = 
"-max and n — n m j n . For the data in figure El we get c « 9.75 x 10 . 
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Numerical simulation - 
Conjectured density - 



Fig. 9. The value distribution for the eigenfunctions of a Seba billiard 



Appendix 

We here show that the distribution of the maximum amplitude Ai(n, L; v) com- 
pletely determines the value distribution of the eigenfunctions on the i th bond, 
which is described by 



1 



JV 



NLi 



/(#>(*)) da 



(51) 



where / is an arbitrary bounded continuous function. Let us in fact consider the 
more general joint distribution, 

1 N [ u 

— — V / F(cosk n (x - Li),v 2 Ai(n,L;v)) dx 
IXL i tT, Jo 



where F is a bounded continuous function in two variables. We obtain the ex- 
pression H51|l for the choice F(t, rj) = f(t^/rj) provided / is even (which, as will 
become clear below, we may assume w.l.o.g.). 

We begin with the special case when F factorizes, i.e., F(t,rf) = fi(t) f^rj) 
where /i, /2 are arbitrary bounded continuous functions. Then 



1 

T 



Li 



i JO 



/i(cosfc„(x — Li)) dx = / /i ( cos(27rx)) dx + 0(fc n 1 ) 



:r ./_ i VI - t 2 



+ 0(fc- 1 ), (52) 



and, by theorem 



i N poo 

/ 2 (tMfn, £;«))-► / f 2 (v)Qv(v)dv 

n—1 ^ 



(53) 
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as N — > oo. Since the mean density of the eigenvalues k n is constant, we have 
E„<jv °( k n X ) = O(logiV) and thus from (52) and (gjj) 



1 N f Li 

lim — — / F(cosfc„(x - L. l ) 1 v 2 A i (n, L;v)) dx 

1 71=1 1/0 



=\LX F ^ )QM 7Bw- (54) 

This holds for functions F — fi f 2 and, by linearity, also for finite linear com- 
binations of such functions. Given any e > we can approximate any bounded 
continuous F from above and below by such finite linear combinations F+ and 
F_ , respectively, such that 

1 f 1 r° r , * , m , N dndt 

Since e can be arbitrarily small, 1|54|) holds in fact for any bounded continuous 
F. 

We can therefore choose F{t,rj) = f(ty/rj) as a test function, and we find 



1 «C30 



with the limiting distribution 



fti(r) = - / Q„(s) 



7T J r 2 y 's — r 2 

The limit u — > cxd can be handled in an analogous way and leads to the 
same formulas for the limit R(r) of R v (r) with Q v (s) replaced by Q(s) in the 
above. It follows from the asymptotic expansion J3J that R(r) also decays with 
an algebraic tail. 
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